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Abstract 

It has recently been shown that there exist stable inhomogeneous neutral 
black strings in higher dimensional gravity. These solutions were motivated 
by the fact that the corresponding homogeneous solutions are unstable. We 
show that there exist new inhomogeneous black string and black p-brane 
solutions even when the corresponding translationally invariant solutions are 
stable. In particular, we show there exist inhomogeneous near-extremal black 
strings and p-branes. Some of these solutions remain inhomogeneous even 
when the size of the compact direction (at infinity) is very small. 
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I. INTRODUCTION 



Higher dimensional generalizations of black holes, such as black strings and black branes, 
have played an important role in string theory. The first solutions of this type were found by 
assuming translational symmetry along the brane jI],H. However it was shown by Gregory 
and Laflamme that many of these solutions are unstable to linearized perturbations 
which are inhomogeneous along the brane. Although it was originally believed that the 
endpoint of this instability (at least in the neutral case) was individual black holes, it was 
recently shown that this is incorrect. Instead, the solution must settle down to stable 
inhomogeneous black strings and p-branes [gj . 

In this paper we show that there exists a large class of new stable inhomogeneous black 
branes. These solutions are unrelated to the Gregory-Laflamme instability and exist even 
when the corresponding homogeneous solution (with the same mass and charge) is stable. 
The new solutions are charged and near their extremal limit. Gubser and Mitra || con- 
jectured that for a black brane with a (noncompact) translational symmetry, the Gregory- 
Laflamme instability exists precisely when the brane is thermodynamically unstable. They 
demonstrated numerically that a certain class of black holes in anti-de Sitter space were 
thermodynamically unstable precisely when they were also unstable against linear pertur- 
bations |7]]. Recently, Reall || provided strong support for this conjecture. In the cases we 
consider, the near extremal homogeneous black branes have positive specific heat, so they 
are expected to be stable. Nevertheless, we find new inhomogeneous solutions which are 
also stable. This shows that there can be more than one stable black brane configuration 
with the same mass and charge. 

We cannot write down the new solutions explicitly. Instead we construct time symmetric 
initial data with an apparent horizon whose area is much larger than that of the homogeneous 
black brane with the same mass and charge. Under evolution, the mass can only decrease 
(through the emission of radiation) and the horizon area can only increase, so the final state 
will still have a horizon with area much larger than the corresponding homogeneous solution. 
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It must be a new stable black brane solution which is not translationally invariant along the 
brane. 

To illustrate our construction of the initial data, consider the simple case of five dimen- 
sional Einstein-Maxwell theory. In addition to charged black holes with S 3 event horizons, 
this theory is known to have translationally invariant charged black strings. To keep the 
total mass and horizon area finite, let us compactify the direction along the string so its 
horizon has topology S 2 x S 1 . Since there is no force between extreme charged black holes, 
one can explicitly write down the solution for a one dimensional periodic array, which is 
equivalent to a black hole in a space with one direction compactified. We can now compare 
the horizon area for the black hole and black string with the same mass, charge, and length 
of the circle at infinity. Interestingly, in the extremal limit, the black string has zero horizon 
area while the area of the black hole remains nonzero. One can now produce new solutions 
by starting with the array of extreme black holes and adding thin neutral black strings to 
connect the horizons. From the compactified viewpoint, the black string wraps around the 
circle and has its ends stuck on the black hole. The total increase in mass can be made very 
small by taking the black string to be very thin, so the configuration is near extremality. 
There is an S 2 x S* 1 event horizon with area comparable to the extremal black hole. This 
configuration is not static and will evolve. But it cannot settle down to the homogeneous 
black string since that has very small horizon area near extremality. Instead, it must settle 
down to a new inhomogeneous near extremal black string. We will show that the homoge- 
neous black string has positive specific heat near extremality, so these new solutions are not 
the result of a Gregory-Laflamme instability. 

The new solutions have at least two surprising properties. First, in constructing the 
periodic array of extreme black holes, the charge Q and size of the circle at infinity L are 
independent parameters. In particular, one can put a large extreme black hole (with horizon 
area proportional to Q 3 ^ 2 ) inside a space with arbitrarily small L. The horizon remains a 
round sphere and is not distorted. Instead, the size of the circle grows as one comes in from 
infinity to accommodate the black hole. Our construction thus gives an inhomogeneous 
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near extremal black string in a space with arbitrarily small L at infinity. It is unlikely 
that a similar statement is true for neutral black strings. (However, see for evidence that 
there exist inhomogeneous neutral black strings even when the corresponding translationally 
invariant solutions are stable.) 

The second surprising property follows from the fact that extremal black holes have 
infinite throats. The thin black string that we add must go all the way down the throat to 
reach the horizon. Thus it has infinite length. This does not contradict our statement that 
the net increase in mass is very small, since the mass is infinitely redshifted near the horizon. 
However the infinite length implies that the area of the black string horizon is infinite. We 
will show that the black string horizon is an apparent horizon for this initial data, so we 
have an apparent horizon with infinite area in a space with finite total energy! Previously, 
the fact that extremal black holes have infinite throats did not seem to have any physical 
consequences since all matter thrown into the black hole reaches the horizon in finite time. 
However, we see that it can have a major effect when one adds the black string. On the 
other hand, the true event horizon will lie outside the apparent horizon and does not have 
to go down the throat. Its area will be finite and probably of order the area of the black 
hole. 

One can do a similar construction for general black p-branes in string theory. We can start 
with a periodic array of extremal black branes if their horizon area is nonzero. Otherwise, 
we can start with a periodic array of near extremal solutions. One can then connect them 
with thin neutral black strings or even neutral black branes. We will see that the horizon 
area can be made larger than the corresponding "smeared" homogeneous solution, so it must 
evolve into a new inhomogeneous near extremal black brane. 

In recent years the Bekenstein-Hawking entropy of near extremal black holes and homo- 



geneous black branes has been understood in string theory. (For a review see ||10|| .) It would 
be interesting to understand the entropy of these new inhomogeneous black branes. We will 
make some preliminary comments about this in section V. 

This paper is organized as follows: In section II, we review and compare the black hole 



and black string solutions of five dimensional Einstein-Maxwell theory. In section III, we 
construct time symmetric initial data describing a periodic array of extreme black holes 
connected by a neutral black string. We also show that this initial data cannot evolve into 
the known translationally invariant solution. In section IV, we extend this construction 
to higher dimensions and argue that there must exist inhomogeneous near extremal black 
p-branes. Finally, section V contains some further discussion, including comments about 
the possible microstates associated with these solutions. 



II. REVIEW OF STATIC SOLUTIONS 



The simplest context to describe the new solutions is five dimensional Einstein-Maxwell 
theory with action 

R 1 
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a,b = 0, ...,4 



(2.1) 



where G is the five dimensional Newton constant, R is scalar curvature, and F is Maxwell 
field. This theory is known to have both electrically charged black holes and translationally 
invariant, electrically charged black strings. The latter can be viewed as resulting from the 
collapse of a cylinder of charged dust. The black holes are described by the five dimensional 
generalization of the Reissner- Nordstrom solution. We are mostly interested in the extremal 



limit which can easily be obtained as follows [11]. Let 



ds' 



-U-\x l )dt 2 + U(x l )Sijdx l dx J 



(2.2) 



Ei = aU^diU (2.3) 

where a = ±(3/167rG) 1//2 . Then the Einstein-Maxwell equations reduce to just the condition 
that U be a harmonic function. If U is the field of a point mass, 
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where [l is a positive constant, then the solution ( |2.2j ) describes an extremal black hole. 
The spatial metric has an infinite "throat" since the proper distance to x % = is stretched 
out infinitely, and near the origin the area of the three-spheres of constant radius is almost 
independent of the radius. In these coordinates, the event horizon is at x % = and has area 

A BR = 2n 2 ^ 2 . (2.5) 

The ADM mass M and charge Q are given by 

m =%* Q=±yf^, (2.6) 

respectively, where we have simply normalized the charge by Q = § EidS 1 . 

Solutions describing several extremal black holes are easily constructed by letting U have 
several point sources. (These are the analog of the familiar Majumdar-Papapetrou solutions 
in 3 + 1 dimensions [|l|,|13|].) To compactify one direction, we let U be the field of a one 
dimensional periodic array of point masses. The resulting metric can be written [11| 

ds\ = -U- 2 {r, z)dt 2 + U(r, z){dr 2 + r 2 dQ 2 + dz 2 ), 
tti x , ku, sinh27rf 

^ ) = 1 + l Lsh2, t -cW < 27) 
where the coordinate z is periodic with period L. The black hole horizon is located at 
r = z = 0, where U diverges. Expanding U near this point yields 

"<'.*> = 1 + ;af? + T& + (2 ' 8) 

So the geometry near the horizon reduces to that of the isolated black hole. For r ^> L, we 
have 

U = 1 + ^ + 0(e- 2 ^/ L ) (2.9) 
rL 

Note that the inhomogeneity in the z direction falls off exponentially for large r. This 
is expected since, from a four dimensional viewpoint, z dependent perturbations act like 
massive fields. The ADM mass M and charge Q of the compactified solution are identical 
with that of the single black hole ( |2.6|) . 
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It is interesting to note that /x and L are independent parameters in this solution: One 
can fit an arbitrarily large charged black hole into a space with one direction compactified 
on an arbitrarily small circle (at infinity). This is possible since the size of the circle depends 
on U . It follows from (|2.9|) that when r ~ L the proper length of the circle is of order /i 1//2 , 



independent of L. 

To obtain the (translationally invariant) charged black string solution we can proceed as 
follows. Consider the following metric and two form 

dsl = e-^^dz 2 + e^^g^dx^dx", /x, v = 0, 3 

F 2 = -F^dx? Adx", (2.10) 
where g^ u is a four- dimensional metric and z is again a compactified extra-dimension with 



period L. As shown in |L4 |, if all functions are independent of z, the action fl2.1|) is reduced 



to the following four-dimensional action 

1 r „ _ ,_ „ I 



S = L fy/=g —-[R - 2 (V0) 2 ] - -e-W^F^ 



d*x. (2.11) 



The electrically charged spherically symmetric black hole solution of this theory is given by 

urn 



ds 2 = -\ 2 dt 2 + X~ 2 dr 2 + fdQ 2 , 

\ 2 = (l- r -±) fi_^ 1/2 

r / V r 



/ r \l/4 

f = r[l--j , (2.12) 

where r + and r_ (< r + ) are free parameters. The Maxwell field F and the "dilaton field" 
are 



and 



F 2 = ± — J-^dtArfr (2.13) 



= 1 - -) , (2.14) 



respectively. Substituting Eqs. (|2.12| ) and ( |2.14j) into Eq. (|2.10|) , we obtain the five- 
dimensional black string solution: 

ds\ = - (l - I±) (l - 1^) dt 2 yY'dr 2 + r 2 (l - ^) dtf + (l - ^'dz 2 . (2.15) 

The event horizon is at r = r + and has 

/ \ 1/2 

A = Anr 2 + L ( 1 - ^— J . (2.16) 

There is a curvature singularity at r = r_. In the extremal limit, r + — > r_, the horizon area 
clearly goes to zero. By a simple coordinate transformation, the extremal solution can be 
put into the form ( [2.2|) , where now U is the field of a line source. The total mass M and 
charge Q is given by 



M=A( 2r+ + r _), Q = ±L J^±^. (2.17) 



AG K T " V G 

In terms of Q and the Hawking temperature 



H Aixr + y r + ^ ^ 

we can calculate the specific heat C defined as 

C(4rH ) = ^§l£). „ 19) 
In the near-extreme limit (r+ — > r_), we obtain the asymptotic form of C: 



C(g,^) = Ji^M-(TH + 0(2f)). (2.20) 



Since the specific heat C is clearly positive in the near-extreme case, the solution (|2.15 ) is 
believed to be stable. 

Equating the mass of the black hole to that of the black string in the extremal limit 



r + = r_, we find fi = Lr + /n. It then follows from Eqs. ( |2.6[ ) and ( 2.17 ) that the charges of 
the two systems are also equal. But we have seen that the horizon areas are very different. 
We will use this fact in the next section to construct initial data for a new inhomogeneous 
near extremal black string. 
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III. INITIAL DATA FOR AN INHOMOGENEOUS BLACK STRING 



In this section we show that one can start with the extremal black hole solution in a 
space with one direction periodically identified, and add a small amount of energy so that 
the horizon changes from having topology S 3 to having topology S 2 x S 1 . The idea is to 
add a thin neutral black string which goes around the circle and has its ends stuck on the 
black hole. This configuration will not be static and will evolve. We do not know the final 
exact static solution, but we will construct time symmetric initial data and then argue that 
it must settle down to a new inhomogeneous near extremal black string. There are in fact 
two different ways to construct the initial data. The first is more explicit, and the second 
is more general but (at the moment) less rigorous. We will concentrate on the first method 
and comment on the second at the end of this section. 



A. Conformal rescaling construction 



Let us consider the following four- dimensional spatial metric: 



ds 4 = u 



= uds 2 , (3.1) 

where u is a function of (r,z,6,(p). z is a periodic coordinate with period L. In the case 
of u — 1, the metric represents the spatial part of the neutral black string solution with 
Schwarzschild radius r = r , while the r = case reproduces a spatial slice in the extremal 
solution ( gg) . 

Assuming time symmetric initial data, the constraint equations reduce to 

R = (3.2) 

= 0, (3.3) 

V is the covariant derivative with respect to the metric ds\. Since the scalar curvature R of 
the metric ds 2 is equal to zero, R can be simply written in terms of u: 



3^, (3.4) 
where V represents the covariant derivative operator of the metric ds 2 . Setting 

Ei = aM _1 VjM, (3.5) 
where a is a constant, the constraint Eq. (|3.3| ) is reduced to the Laplace equation 

V 2 u = 0. (3.6) 
Substituting Eq. (|3.5|) into Eq. ( ^2|) , we can obtain a as 



a = ±V isg- (3 - 7) 

So Eq. (|3.6|) is the only constraint equation for the time symmetric initial data. 

Let u be the solution with a point source at 9 = 0, r = r , z = 0. This corresponds 
to adding an extreme black hole at the north pole of the black string. Introducing a new 



coordinate r* = J r r o dr/Jl — r /r, we can make an analytic continuation to the region r* < 0. 
Hence this initial data has a symmetry r* «-> — r*. Since the initial data is time symmetric, 
the outgoing null expansion 8 + for any closed three-dimensional hypersurface is given by the 
extrinsic curvature of the hypersurface on the initial data. This means that r* = (r = ro) 
is an apparent horizon. To estimate the area, let us consider the local geometry around the 
north pole located at r* = and 9 = 0. Since the spatial metric ds 2 is locally flat, ds 2 is 
approximately 

ds 2 ~ dx 2 + df + dz 2 + drl, (3.8) 

where (x, y) are local orthogonal coordinates around the north pole. The approximate 
solution of u becomes 

u ~ for R 2 = x 2 + y 2 + z 2 + r 2 < r , (3.9) 

i? 2 

where /x is a positive constant proportional to the charge. The area of the apparent horizon 
around the north pole is approximately given by 
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A ~ 



dR 



(3.10) 



oo. 



R 



r*=0 



This shows that the apparent horizon (r = tq) has an infinite area! Since the apparent 



At first glance, the area of the event horizon also seems to be infinite because the event 
horizon is outside the apparent horizon. This speculation is obviously wrong because outside 
of the apparent horizon u is finite everywhere on the hypersurface. Although the exact 
location of the event horizon is unknown, we can still produce a lower bound on its area. 
The idea is to consider a region of the space and show that all surfaces passing through this 
region have area larger than some constant. Our estimate will be crude, but sufficient to 
establish the existence of new inhomogeneous solutions. 

Adding the extreme black hole at the north pole of the black string horizon produces 
initial data which is approximately spherically symmetric for r ^> r . It will be convenient 
to have initial data which is approximately spherically symmetric even for r ~ 2tq. This 
can be obtained by distributing a large number N of extreme black holes each with charge 
2Q/N > uniformly on the sphere (r = tq, z = 0) 0. We now construct a region such that 
the effects of both the compactified extra dimension and the tiny black string are negligible. 
To this end, let = Nr . By choosing r small enough we can arrange for < L/N. 
Introducing a new radial coordinate R 2 = r 2 + z 2 and angular coordinate x such that r = 



1 Strictly speaking, the topology of the apparent horizon is S 1 x S 2 minus a point p. However, 
one can modify the horizon so that the topology becomes S* 1 x S 2 by covering the neighborhood 
of p with a very shallow cap. As easily checked, the outgoing null expansion on the cap cannot be 
positive since u rapidly decreases toward the outside. This produces a compact apparent horizon 
with arbitrarily large (but finite) area. 

2 Because our initial data has two asymptotically flat regions, half of the total electric flux reaches 
each one. Therefore, we put 2Q on the sphere (r = ro,z = 0) as a total charge so that the total 
flux for one asymptotically flat region is just Q. 



horizon has S 1 x S* 2 



topology, the event horizon also has S 1 x S 2 topology []. 
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.Rcosx, z = i?sinx, define region B by (R,£) G [V%r b < R < L/N, —n/A < x < 7r/4] (see 
Fig. [I]). Then, the approximate solution in Eq. (|2.8| ) is valid everywhere in B. Also define 
region A by the almost all spherically symmetric region enclosed by the u = const, surface 
intersecting a point (r = r b , z = r b ) (see Fig. Since the electric field normal to each 
u = const, surface should be positive, i.e. Eiri 1 ~ — (Viu)n % > 0, (where n l is a unit 
outward normal vector for each u = const, surface), u must be increasing as one moves in. 
So the following lower bound on u is obtained: 

U\A>U\(r,z)=(r b ,r b )=U b 3 I 3 - 11 ) 

r b 

by Eq. Q (see also Fig. |1|). 

Now, one can estimate a lower bound Aq for the area of the event horizon. Fig. [I] shows 
typical cases for the event horizon passing through our regions A or B. In general it is 
possible that the horizon is more complicated and wavy. In that case the area would be 
greater than the lower bound of the typical cases. So, it is enough for the estimation of A 
to consider only the typical cases. Recalling the spatial metric Eq. ( |3.i|) , let us first consider 
the case that the event horizon intersects the region A. Then, the area of the event horizon 
in this region A is roughly 

,2 ir >/2 

where we used r /r b = 1/N. Note that this lower bound is independent of r . In the case 
that the event horizon intersects the region B, we can easily obtain a lower bound on Aq. 
Eq. (|2.8| ) is the approximate solution for the region B, and we will assume the constant term 
is negligible. (If necessary, we can choose iV so that fi/L 2 ^> 1/N 2 to insure this.) Then 
the area of each u = const, surface is almost /i 3//2 . So, if the orbit of the event horizon is 

R = R(x), 

A ~ 4vr J~[(uR 2 f> (l + (^j j cos 2 X d X 

7T 

> 4vr f \ {uR 2 ) i cos 2 x dx ~ ^ , (3.13) 
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A ~ ^(v^r) 2 ^^) > [j A {^T b r) 2 ^r b dz^ > ^ 2 = (3.12) 



where a dash denotes the derivative with respect to x- We used Eq. ( |2.8|) to derive the final 
value. Finally, suppose the event horizon passes outside of the region B, and consider the 
subset with —L/N < z < L/N. In this region u is given by ( |2.7| ). From this explicit solution 
one can show that ur 2 > 11. So the proper area of the two-spheres in this region is at least 
Li. Since the length in the z direction is at least L/N (since u > 1 everywhere) we obtain 
the lower bound Aq > 11L/N . Thus, we can get the lower bound for the area of the event 
horizon as 

A o>min|^, / ,-|, (3.14) 

which is independent of r . 
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FIGURES 




z=0 



r 

FIG. 1. Constant u surfaces and possible cases for the event horizon are depicted on the (r, z) 
surface. Each point represents a two-sphere. For the region B, the event horizon would be more 
convex than the constant u surfaces because the outgoing null expansion on those surfaces is almost 
zero due to the infinite throat. (The outgoing null expansion along the event horizon should be 
strictly positive when the spacetime is not static). 

B. Evolution of initial data 

The initial data constructed above must evolve since, e.g., the surface gravity is not 
constant over the horizon. It is large near the thin black string and small near the extremal 
black hole. More physically, one would expect the extremal black hole to start swallowing 
up the string. What will it evolve to? We will assume cosmic censorship holds, since a 
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generic violation of cosmic censorship in five dimensions would be an even more surprising 
conclusion than we find below. We will also use the result of || that cycles on an event 
horizon cannot shrink to zero size. This means that if the event horizon is initially S 2 x S 1 , 
it must remain S 2 x S 1 . It cannot pinch off. 

The ADM mass of our initial data set is easily computed from (|3.1|). Since u is a solution 



to Laplace's equation in the black string background with a point source, the asymptotic 
form of u can be directly related to the source by integrating the equation over all space. This 
means that the asymptotic form is independent of tq and the same as (|2.9| ): u ~ 1 + irfi/Lr. 
The ADM mass is thus 

M ADM = M+ r ^. (3.15) 

where M is related to /i by ( p.6|) . The total charge Q is also related to ji by Q2.6Q . So the 
mass is above the extremal limit just by the second term which can be made very small 
by choosing r small. Under evolution, the mass can only decrease since energy can be 
radiated away to infinity. Since the charge is conserved, the final state must be even closer 
to extremality. On the other hand, the horizon area can only increase during the evolution 
by the usual area theorem. Thus the final state must be a black string with horizon area 
larger than the initial fl3Tj ). But a homo geneous black string (|2.15|) with small excess 
energy tqL/2G above extremality would have r + — r_ = 4r . Thus its horizon area (|2.16|) 



1/2 

would be A oc r . By choosing tq sufficiently small, this can clearly be made much smaller 
than the initial horizon area. So the initial data cannot settle down to the homogeneous 
black string. It must approach a new stable, inhomogeneous near extremal black string. 

We have been discussing classical solutions, but if one has to take tq less than the Planck 
scale (or string scale) the solutions would not be physically interesting. Quantum effects 
would certainly be important near the horizon. Fortunately, it is easy to see that this is 
not necessary. For our lower bound on the horizon area, iV ~ 100 should be more than 
adequate. So as long as /z and L are much larger than the Planck scale, we can choose tq 
greater than the Planck scale and still deduce the existence of new solutions. 
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It is plausible that the final configuration will resemble a large near extremal black hole 
with a near extremal black string going around the circle^. This configuration would have 
low surface gravity everywhere. If one starts with this inhomogeneous black string and takes 
its extremal limit, it is plausible that the solution degenerates to an extremal black hole with 
a (singular) extremal black string going around the circle. This solution can be constructed 
explicitly. Recall the general form of static extremal solutions (|2.2|) . The solution we want 
is obtained by letting U be the harmonic function with a line source and a point source 
superposed on the line. The line source produces the extremal black string, while the point 
source reproduces the extremal black hole. In fact, there is a large class of extremal solutions 
of this type, since one can specify the line density arbitrarily. It is not clear whether there are 
near extremal analogs of all of these extremal solutions, or just a few. If all of these extremal 
solutions can be obtained as the extremal limit of some black string, then the structure of 
black string solutions would be very rich indeed. It would depend on an arbitrary function. 

C. Gluing construction 

The previous construction of time symmetric initial data is very explicit, but restricted to 
situations where one wants to add black strings to extremal black holes. In the next section 
we will want to add black strings to near extremal black p-branes. To do this we need a 
more general construction which we describe here. This construction is less explicit (and, 
at the moment, less rigorous) but can be used to insert a thin black string into essentially 
any initial data set. 

Intuitively, it seems clear that one should be able to take any smooth initial data set 
and insert a very small Schwarzschild black hole. This is because every curved space is 
locally flat, and the Schwarzschild black hole is approximately flat after a few Schwarzschild 
radii from the black hole. (Near the black hole there is large curvature, but the constraint 



3 We thank S. Ross for suggesting this. 



16 



equations are satisfied.) This has indeed been proven in 3 + 1 dimensions in a recent paper 
(T7|j . We would like an analogous result for black strings in 4 + 1 dimensions. In other words, 
it should be possible to remove a neighborhood of a geodesic in an initial data set and glue 



in a thin black string. This possibility is supported by the fact that the argument in (17 
can be easily generalized to show that one can glue small black holes in 4 + 1 dimensions 
| HI . But if one adds a sequence of closely spaced black holes, there will be a single apparent 



horizon enclosing all of them which is essentially a black string. 

In the case of interest here, we start with the exact solution describing a large five 
dimensional extremal black hole with the z direction compactified. The curvature on a 
static slice is everywhere small. The z axis comes up the infinite throat of the extreme black 
hole, goes around the circle and then back down the infinite throat. We now want to remove 
a neighborhood of this curve and glue in a thin black string, i.e., 4D Schwarzschild cross 
R. The fact that the initial data contains an electric field should not be a problem. Let the 
electric field along the z axis be E z (z). Since the electric field is parallel to the z axis, when 
we glue in the black string, we can add an electric field (pointing in the R direction) equal 
to E z (z). 

Although we cannot write down this initial data explicitly, we can write down a metric 
which describes the solution asymptotically down the throat of the extreme black hole. 
Before we add the black string, the asymptotic solution is just 5 3 xR with a constant electric 
field directed up the throat. So the constraint equation is just that the scalar curvature is 
a positive constant. Adding the black string on either side of the throat, corresponds to 
adding two small black holes at the poles of the large S 3 . But this is exactly the geometry of 
a static slice of Schwarzschild de Sitter in 3 + 1 dimensions. This metric has constant positive 
scalar curvature by virtue of the Einstein constraint equation with a positive cosmological 
constant. So the initial data down the throat of the extremal black hole with a black string 
added is just the product of a line and a static slice of 3 + 1 Schwarzschild de Sitter. The 
electric field is still constant and directed along the line. 

The resulting initial data differs from the one obtained by the conformal approach in the 
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following respect. In the conformal approach, the size of the black string apparent horizon 
is rescaled by u. Thus, near the extremal black hole it becomes much larger than r . In 
fact, it is easy to see that the r = surface intersects the extremal black hole horizon at 
the equator. So the black string has grown to a size given by Ml In the gluing approach 
this does not happen. The size of the black string stays small, O(r ), all the way down the 
throat. 

IV. GENERALIZATION TO P-BRANES 

A. Starting with extremal branes 

The above construction can be applied to a variety of higher dimensional situations to de- 
duce the existence of inhomogeneous near extremal p-branes. First we place the example we 
have already discussed in the context of string theory. Five dimensional Einstein-Maxwell 
theory arises from ten dimensional (Type II) string theory by compactifying on T 5 and 
considering three fundamental charges. The standard choice is Dl, D5 and momentum^. 
However, due to the momentum, the ten dimensional metric has an ergoregion which com- 
plicates the construction. This can be avoided by applying a T-duality, so that the solution 
contains DO, D4, and fundamental string charge. If we let the charges all be different, the 
extremal solution has the following (string frame) metric [f2T| : 

ds 2 = - H 1/2 Hl 1/2 H{ l dt 2 + H^H^idyidy') (4.1) 
+ Hl /2 Hl /2 (H{ l dx 2 + dr 2 + r 2 dn + dz 2 ) 

where 



4 Starting with this extremal solution, one can add a traveling wave along the Dl direction fl^| . 
However the local horizon geometry remains homogeneous and is independent of the wave profile 
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are three harmonic functions, and yi are coordinates on T 4 . The 4-branes are wrapped 
around this T 4 , and the fundamental strings are wrapped in the x direction which we also 
assume is compact. The branes are all distributed uniformly over these compact directions. 
There are four directions transverse to all the branes labelled by z, r, 8, 0. There is an event 
horizon at the origin of these coordinates, z = r = 0, with nonzero area. Since this event 
horizon is extended in the x, t/i directions, the solution describes an extremal black five- 
brane. Its horizon has topology S 3 x T 5 . If we set yUo = /^i = A*4? then the solution ( f4.1| ) 
reduces to the product of (|2.2|) and a flat T 5 . 

Since Hi are just harmonic functions, we can compactify the z direction by taking a 
periodic array as before. The homogeneous solution with the same mass and charges can be 
obtained by choosing the harmonic functions to be independent of z, Hi — 1 + fx.Jr. The 
near extremal solution is obtained by adding a factor of (1 — r /r) to g tt and (1 — r /r) _1 to 
g rr as usual. It is easy to verify that the near extremal homogeneous solution has positive 
specific heat, so by the Gubser-Mitra conjecture it should be stable. Nevertheless, we now 
show that there exist near extremal solutions which are inhomogeneous in the z direction. 

To begin, note that the horizon of the homogeneous solution becomes singular (zero 
area) in the extremal limit. This is expected since if the solution is homogeneous in the z 
direction also, it can be dimensionally reduced to a four dimensional black hole. But it is 
well known that extremal black holes in four dimensions will have nonzero area only if there 



are four nonzero fundamental charges p2],|23[. Since we have only three nonzero charges, the 
four dimensional extremal black hole must have zero area. 

Now, starting with the array of extremal black five branes in ten dimensions, one has 
several possibilities to construct initial data that will evolve into a new inhomogeneous black 
brane solution. The basic idea is the same as before. By adding a small amount of energy 
we can cause the horizon to extend in the z direction. Since the initial horizon area is much 
larger than the corresponding near extremal homogeneous solution, the final state must be 
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something new. 

In ten dimensions, there are neutral black p-branes which are just the product of R p 
and the (10 — p)-dimensional Schwarzschild solution. We can use any of the solutions with 
1 < p < 6 to connect the extremal black five-branes in the periodic array. For example, 
one possibility is to add a thin neutral black string in ten dimensions which goes around 
the z direction and ends on the black five-brane. This will produce an apparent horizon 
with topology S 3 x T 5 with aS 7 xR handle addedQ. Under evolution, the topology of the 
apparent horizon can change, so the final topology of the event horizon may be different. But 
the solution cannot settle down to the known homogeneous solution. Another possibility is 
to use a thin neutral black six brane to connect the extremal elements in the array. (This 
is just the lift of the black string we added in five dimensions in the previous section.) This 
will produce an apparent horizon with topology S 2 x T 6 , which is the same topology as the 
horizon in the near extremal homogeneous solution. But once again, the initial area is too 
large for the solution to become homogeneous. Clearly, one can do a similar construction 
with any of the thin neutral black p-branes with 1 < p < 6. 

One should be able to use either the gluing approach or a generalization of the conformal 
approach (where one rescales different parts of the spatial metric by different conformal 
factors) to construct suitable initial data. 

B. Starting with near-extremal branes 

The above solutions contained three nonzero charges in ten dimensions. If one starts 
with just a single charge, then the horizon area always goes to zero in the extremal limit. 
For example, consider a near-extremal D3-brane wrapped on a T 3 of volume V: 



5 That is, remove two small balls from S 3 x T 5 . The boundary of each ball is an S 7 . Now identify 
these two boundaries. 
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ds 2 



H. 



-1/2 



(P) 



dt 2 + dyidy 1 



+ H 1 3 /2 (p) 



1 d P 2 + p 2 dn 5 



(4.3) 



where Hs(p) = 1 + Pzl P ■ The horizon is at p = p and has area A = h 3 PqVH 3 (p, 



vl/2 



This area goes to zero in the extremal limit (po = 0, fixed) so if we start with a one 
dimensional array of extremal three-branes and then add neutral black strings, the area of 
the apparent horizon will still be small and not obviously larger than the corresponding 
homogeneous solution. However, we can still show the existence of inhomogeneous solutions 
if we start with an array of near-extremal 3-branes. Such solutions are not known explicitly, 
but if the separation between the branes L is large compared to po> it is plausible that one 
can approximate the solution by starting with ( }4.3|) and letting H be the harmonic function 
on flat R 6 with a one dimensional array of point sources separated by L. Near each brane, 
the metric will reduce to the single brane solution (|4.3|) , and for r = 0, the metric reduces 
to the known array of extremal three branes. Let us label the direction along the array by 



z. 



We now want to compare the horizon area of this near extremal array with that of a 
solution with the same mass and charge which is uniformly smeared over the z direction. 
That solution is easily obtained by starting with the near extremal D4-brane and applying 
T-duality along one direction of the brane: 



ds 2 



Hs 1/2 (r) 



l--^\dt 2 + dyidy 1 



dr 2 + r 2 dVL A + dz 2 



(4.4) 



where H 3 (r) = 1 + pz/r 2 ". One can easily check that this solution has positive specific heat 



and should be stable to small perturbations. It turns out p4j that the array has larger 
horizon area than the homogeneous solution whenever L is larger than r . (Unlike the 
previous cases where we could let L be arbitrarily small, here we require L > r .) We are 
now in a similar situation as before and can deduce the existence of new inhomogeneous 
solutions. The simplest possibility is to connect the elements of the array with a neutral 
black 4-brane, i.e. the product of 6D Schwarzschild and T 4 . The 4-brane wraps the same T 3 
as the 3-brane, plus the z direction. Since the curvature of the near extremal 3-brane is small 
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near the horizon if the charge is large, it should be possible to construct this initial data 
using the gluing construction. By exactly the same arguments as in section 3, this initial 
data will evolve, but it cannot settle down to the homogeneous smeared 3-brane solution 
since that solution has a horizon which is too small. It must settle down to a new smeared 
black 3-brane, which is near extremal but inhomogeneous in the direction perpendicular to 
the branes. 

Another possibility is to add a ten dimensional neutral black string to connect the ele- 
ments of the array. The apparent horizon would now have the topology of T 3 x S 5 with a 
S 7 x R handle attached. As before, the final topology of the event horizon could be different, 
but it could not settle down to the homogeneous solution. One could also glue in neutral 
black 2-branes or 3-branes. 

We have started with large near extremal 3-branes. But one could equally well start with 
any near extremal brane such that the curvature is small near the horizon. This includes 
the M2 and M5 branes in eleven dimensions. 



V. DISCUSSION 

We have shown that there exist a large class of inhomogeneous extended black holes 
in higher dimension. Unlike the solutions discussed previously 0, these solutions exist 
even when the translationally invariant black branes are stable. Some of these solutions 
remain inhomogeneous even when the length of the compact direction at infinity goes to 
zero. As one consequence, this means that for certain values of the mass and charge, there 
is more than one stable solution. In other words, the four dimensional black hole uniqueness 
theorems do not extend to higher dimensions with horizon topology S p x T q . It should be 
noted that even for the vacuum Einstein equation, it has recently been shown that there is 
no uniqueness theorem. In five dimensions, in addition to the rotating black hole, there is 
another stationary solution with the same mass and angular momentum which describes a 
black string in the shape of a large ring PS[. The ring is rotating and the centrifugal force 
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balances the gravitational attraction. 

The inhomogeneous solutions we have discussed have higher Bekenstein-Hawking en- 
tropy than their homogeneous counterparts. So one might expect that even though the 
homogeneous solution is classically stable, it might quantum mechanically tunnel to the 
inhomogeneous state. We cannot yet estimate the probability for this transition, but expect 
it to be very small for any macroscopic solution. 

An important open question is whether one can give a microscopic description of the 
entropy for these new solutions. For many of the homogeneous near-extremal solutions such 
a description has been obtained in string theory by taking the limit as the string coupling 



g goes to zero ||10|| . The mass and charge remain but no longer gravitate. So the spacetime 
becomes flat. The charge is represented by D-branes and the excess energy is described by 
fields on the brane. The number of states turns out to be the exponential of Sbh — A/ 4 
where A is the area (in Planck units) of the horizon which appears when one increases g. Is 
there a similar description for the inhomogeneous solutions? 

The answer is not yet clear. One obvious problem is that we do not yet have the static 
inhomogeneous solution explicitly to compute its entropy. However, given the intuition about 
what the final state looks like, we can begin to construct a microscopic model. Consider the 
case of near extremal 3-branes connected by a neutral thin 4-brane. We expect that it will 
settle down to a near extremal 3-brane connected by a smeared 3-brane. At weak coupling, 
this corresponds to a large number N of 3-branes all localized at the same point on the circle, 
together with a smaller number N' of 3-branes uniformly spread around the circle. The states 
of the localized 3-branes are described, as usual, by a 3 + 1 dimensional U (N) gauge theory. 
The states of the smeared 3-brane are most easily counted in weak coupling by applying 
T-duality ||24j| . Then it becomes a near extremal 4-brane with nontrivial flat connection, and 



the low energy excitations are described by a U(N') gauge theory in 4 + 1 dimensions. One 
might hope that the entropy could be reproduced by combining both systems. Of course, 
at this level, one does not see any sign of inhomogeneity. That must arise from interactions 
between these states. One can view the existence of these new supergravity solutions as 
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giving a prediction for what happens at strong coupling. 
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